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Abstract. Non-colliding Brownian particles in one dimension is studied. Brownian particles 
start from the origin at time and then they do not collide with each other until finite time T. 
We derive the determinantal expressions for the multitime correlation functions using the self-dual 
quaternion matrices. We consider the scaling limit of the infinite particles N ^ oo and the infinite 
time interval T — > oo. Depending on the scaling, two limit theorems are proved for the multitime 
correlation functions, which may define temporally inhomogeneous infinite particle systems. 



1 Introduction 

We consider the process X{t), which represents the system of Brownian motions in one 
dimension all started from the origin and conditioned never to collide with each other up 
to time T. If we take the limit T — > oo, the system becomes a temporally homogeneous 
diffusion process Y{t), which is the Doob /i-transform |3] of the absorbing Brownian motion 
in a Weyl chamber 

= s X = (xi, ^2, . . . , xn) G R^; Xi < X2 < ■ ■ ■ < X]sf 



with harmonic function /;,Ar(x) = Y[i<i<j<Ni^j ~ ^i) |H]- By virtue of the Karlin-McGregor 
formula [121 CSI, its transition density /Ar(t, x, y) from the state x to y in R^ in time period 
t > is given by 

/7v(t,x,y) = ^ det (ptixi,yj)), 



where pt{x,y) = 7=6 '^^ y)^/^*, Qn the other hand, if the non-colliding time interval T 



/2nt 



remains finite, the process X(t), < t < T, is temporally inhomogeneous [T^ . 

We notice an integral formula found in Harish- Chandra jH], Itzykson and Zuber ^UI, and 
Mehta HE!, 



= c aU exp 



-^tr(X - U^YUy 



^det (pj(a;i,|/j)) 

/lAr(x)/lAr(y) 

with c-^ = (27r)^/2^^V2]^^^+ir(i), where X and Y are the X x iV diagonal matrices, 
Xij = Xi6ij,Yij = yi6ij, and the integral is taken over the group of unitary matrix U of 
size N. This equality implies that the non-colliding Brownian motions such as X{t) and 
Y{t) can be described by using the eigenvalue-statistics of Hermitian random matrices in 
Gaussian ensembles [ISj. In earlier papers^! E], it was shown that Y{t) is identified with 
Dyson's Brownian motion model with /3 = 2 and the particle distribution is expressed by 
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the probability density of eigenvalues of random matrices in the Gaussian unitary ensemble 
(GUE) with variance t, while y j^^^^it) coincides with the distribution of eigenvalues of 



random matrices in the Pandey-Mehta ensemble fHl 123 with a = y-^j^, and this tempo- 
rally inhomogeneous process exhibits a transition from the GUE statistics to the Gaussian 
orthogonal ensemble (GOE) statistics as the time t goes on from to T. 

It is known that the eigenvalue distributions of Hermitian random matrices have deter- 
minantal expressions. For instance, in the GUE, the probability density of eigenvalues is 
expressed by 

Pn{xi,X2, . . . ,xn) = ^ det {KN{xi,Xj)), 

iV ! i<i,j<N 

with K^ix^y) = YIi^q Vii{x)ipi>{y), where 

(1.1) = ^e-^'"H,{x) 

with the £-th Hermite polynomial Hi{x) and hi = ^Jiil^ i\. By the orthogonality of ^kix)-, 
we can prove the equality 

(1.2) / det Kt^ixi.x^dxf^: = [N - N' ^V) det K^^ixi^x.) 

J ^<i,j<N' ^ l<i,j<N'-l ■' 

for any 1 < A^' < A^. Such integral property enables us not only to obtain determinantal 
expressions for correlation functions, but also to argue the N —>■ oo limit of the system 
by studying the large asymptotic of the function Kiy{x,y). With proper scaling limit, 
determinantal point processes with sine-kernel and Airy-kernel are derived. See j2ZI and 
references therein. 

In the present paper, we derive the determinantal expressions of the multitime correlation 
functions for the process X{t). Our aim is to prove limit theorems of the multitime correlation 
functions in the scaling limits of infinite particles N oo and infinite time interval T ^ oo. 
Depending on the scaling, we derive two kinds of limit theorems, one of which provides a 
spatially homogeneous but temporally inhomogeneous infinite particle system (Theorem 1), 
and other of which does the system with inhomogeneity both in space and time (Theorem 2). 
We remark that it is easier to prove the limit theorems for Dyson's Brownian motion model 
Y{t). Corresponding to Theorem 1, we will obtain the multitime correlation functions of the 
homogeneous system, which coincides with the system studied by Spohn [2H], Osada 
and Nagao and Forrester j2l]. Similarly, corresponding to Theorem 2, an infinite system 
with spatial inhomogeneity will be derived, which is related with the Airy process recently 
studied by Prahofer and Spohn ^22j and Johansson P3]. 

One of the key points of our arguments is that, in order to give the determinantal ex- 
pressions for the correlation functions for the present processes, we shall prepare matrices 
with the elements, which are neither real nor complex numbers, but quaternions 

q = qo + qiei + ^262 + gscs e Q 
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with gj G C, < i < 3, in which the four basic units {1, Ci, 62, 63} have the following 2x2 
matrix representations, C : Q 1— > Mat2(C); 




The dual of a quaternion q is defined by = go — J2i.=i Qi^i^ ^ quaternion matrix 

Q = lij ^ Q5 its dual matrix = ((Q^)ij) is defined to have the elements {Q^)ij = gjj. 
Following Dyson's definition of the quaternion determinant for self-dual matrices [51 IT71 ITH]. 
we can give the quaternion determinantal expressions having the similar properties to ()1.2|) 
for arbitrary multitime correlation functions for X{t) (Theorem 3). As briefly reported in 
|2n!, the present results can be regarded as simple applications of the results given in Nagao 
and Forrester [22] and Nagao [20] for multimatrix models, and in Forrester, Nagao and 
Honner jB] for the asymptotic of quaternion determinantal systems, here we give, however, a 
self-contained explanation for all the formulae and calculus developed in the random matrix 
theory, which are used to prove our limit theorems. 

The theorems proved here mean the convergence of processes in the sense of finite dimen- 
sional distributions. As argued in Prahofer and Spohn [^ and in Johansson ^1], tightness 
in time should be confirmed. 



2 Statement of Results 

For a given T > 0, we define 
(2.1) ^^(s,x;t,y) 



Tr ^ -.N /iv(^--5,x,y)A/'7v(T-t,y) 



for < s < t < T, X, y G R<, where A/Ar(t, x) = Jj^jv fN(t,yi,'y)dy, which is the probability 

that a Brownian motion started at x G does not hit the boundary of up to time 
t > 0. The function (7^(5, x; t, y) can be regarded as the transition probability density 
from the state x G R^ at time s to the state y G R^ at time t, and associated with the 
temporally inhomogeneous diffusion process, which is the N Brownian motions conditioned 
not to collide with each other in a time interval [0,T]. In jl4[ [T3] it was shown that as 
|x| — i> 0, 5(^(0, X, t, y) converges to 

N 

(2.2) ^^(0, 0, t, y) = C{N, T, t)hM^N{T - t, y) l[pt{0, Vi), 

i=l 

where C{N,T,t) = vr^/^ (jlf^^r{j /2)^ ^ j.N{N-i)/i^-N{N-i)/2 _ ^^len the diffusion process 
X{t) starting from can be constructed. 
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We denote by X the space of countable subset ^ of R satisfying tl(^ H fT) < oo for any 
compact subset K. We introduce the map 7 from f^" defined by 7(0:1, 0:2, . . . , a;„) = 
{xj}"^]^. Then = 7X(t) is the diffusion process on the set X with transition density 

function 'g^{s,^; t,!]), < s < t < T: 

gjf{s, x; t, y), if s > 0, jj^ = [jr/ = A^, 
g'Nis,C;t,v)= { 9%{0,0;t,y), if s = 0, ^ = {0}, ti^ = iV, 



0, 



otherwise, 



where x and y are the elements of with C, = 7X, rj = 7y. For x^^ G R^, 1 < m < M + 1, 



(m) 



and A^' = 1, 2, . . . , A^, we put x^'' = (^x'(^\x^P\ . . . , xj^^ j and = 'j^)^'' ■ For a given 

time interval [0,T], we consider the M intermediate times < ti < ^2 < ■ • • < < 7"- 
Then the multitime transition density function of the process is given by 



M 



(2.3) 



ml tm+ly Cm+l) ■> 



m=0 



where, for convenience, we set to = 0, tM+i = T and C,^ = {0}. From ()2.1|1 and ()2.2|) we 
have 

(2.4) pl{h,^^;t,,^^; . . . ;W,eM+i) = C{N,T,t,)h^ (x«) sgn (/i^ (x^^^^^^)) 

N M 



(m) (m+1) 



i=l 



m=l 



For a sequence {Nm}rn=i of positive integers less than or equal to A^, we define the 
{Ni, N2, . . . , A^Af+i)-niultitime correlation function by 



(2.5) (^ti, ^i^;t2, ^2^; . . . ; ^m+i, ^a/+i^ 



M+l 

n 



TV 



lim— 1 m=l 



n c;a:^V?;(ti,ef;t2,e2^... 



iV ^ 



i=Nm+l 



We will study limit theorems of the correlation functions as A^ — 00. First, we 
consider the case T/v = 2A^. Let 



{s,X]t,y) 



r 1 .1 



TT 



d\ cos(A(a;-y))e-^'(*-^)/2 

sin(x — y) 
Ti{x - y) 



if s > t 



if s = t 



1 r°° 



d\ cos(A(x - i/))e-^'(*-^)/2 if s < ^ 



[s,x;t,y) 



1 

/ (iA A sin(A(x — ?/))e 
Jo 



(s+i)A2/2 



~ 1 1 

x; t, y) = -z I dXj sin(A(x - y))e(^+*)^'/^ 



TT 



And let q'^'"'(x, ?/) be the quaternion, whose 2x2 matrix expression is given by 



^{smj X', Sn, y) I(sm, a^; Sn, 
3{sm,x;sn,y) E>{sn,y;sm,x) 



^(2) 



(Af+1) 



Let M > 1 and {Nm}„:=l be a sequence of positive integers. We denote by Q ( x^^, x^^, . . . , js.^^^^^ 

the self-dual X]m=i x X]m=i quaternion matrix whose elements are q™'" 
l<i< Nm, I <j < Nn, I <m,n < M + 1, that is, 



.(1) ^(2) M+1 



5,A/+l,Af+l /"^(Af+l) ^(M+1) 



with blocks of x A^n quaternion matrices 



,n I (m) (n) 



l<i<N^,l<j<Nr, 



for l<m, n<M + l. 

For an NxN self-dual quaternion matrix Q, the quaternion determinant TdetQ is defined 
by Dyson [5j as 

e{n) 

Tdetg = J2 i-lf-'^''^ n ^-^^^'^ ■ ■ ■ ^'^^ 

vtESat 1 

where £(7r) denotes the number of exclusive cycles of the form {a^b—>-c—>----^d—>-a) 
included in a permutation n G Sn- 

Theorem 1. Let Tj^ = 2N. For any M > 1, any sequence {A^m}m=i of positive integers, 



and any strictly increasing sequence {sm}m=i of nonpositive numbers with sm 



0, 



hm pI- (t^ + si, e^; + ^2""^ • • • ; ^iv, eSt 



TV 

Tdet 



(1) (2) ^(A/+l) 



Remark 1. The above system is spatially homogeneous, since all elements of the quaternion 
determinant are functions of difference of positions, x,-™'^ — x^"''. This expresses the bulk 
property of our infinite particle system. When M = 1, the present system is equivalent with 
the N ^ oo limit of the two-matrix model reported by Pandey and Mehta ES] • In the 
system defined by Theorem 1, if we take the further limit such that Sm —oo with the time 
difference s„ — Sm fixed, 1 < m,n < M, then D(sm, x; s„, y) oo, I{sm, x; s„, y) — >• 0, while 
the product D(sm, x; Sn, y)I{sm, x; Sn,y) — > 0. Therefore, we may replace D and I by zeros in 
this limit, and the quaternion determinant Tdet Q (x^j,x^^, . . . ,xlt^M will be reduced to 



^^2, . . . , y^N2) ^^^^ elements a""'"" (^x'f^\ x'f''^ = 

Hence, we obtain a temporally and spatially homogeneous system, 
whose correlation functions are given by 



an ordinary determinant det A ( xl!,'', xl^^ 



(m) (n) 



p' ^2,^2'"^ • • • ; SM,^Z'') = det A (x(i;,x(') 



MO 



Such a homogeneous system was studied by Spohn j2H] , Osada [24 and Nagao and Forrester 
j2Tj as an infinite particle limit of Dyson's Brownian motion modeH4,. 

Next, we consider the case that T/v = 2N^^^. In order to state the result, we have to 
introduce the following functions. Let Ai(z) be the Airy function: 



(2.6) 

For s,t < and x,y eH, we put 
V{s,x;t,y) -- 



ZTT 



X{s,x;t,y) 



POO J 

/ dX e^^/2Ai(a; + A)-^ {e'^'''M{y + A)} 
Jo "-^ 

poo 1 

- / d\ e*^/2Ai(y + A)-^ {e'^^'^Mix + A)} 

f^OO POO 

/ dXe'^/^Ai{y-X) / dX' e'^' Ai{x - X' 
Jo J\ 



and 
with 



S{s,x;t,y) 
V{s,x;t,y) 



dX e'^'/^Aiix - X) j dX' e'^''^ki{y-X') 
S{s, x; t, y) = S{s, x; t, y) - V{s, x; t, y)l{s < t), 



1 

dX e^'-'^^/^Ai{x + X)Ai{y + A) + -Ai{y) / dX e'^/^Ai{x - A), 



dX e(*-")^/2Ai(x + A)Ai(y + A), 
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where l(s < t) = 1 if s < t, and = otherwise. And let q™''"(a;, y) be the quaternion, whose 
2x2 matrix expression is given by 



C(q™'"(x,i/)) 



Si^Sfni y) X, Sfii y) 



Let M > 1 and {A^^j^+i be a sequence of positive integers. We denote by Q I x 

the self-dual J2m=i ^ X]m=i quaternion matrix whose elements are q'"'" (^x, 
l<i< Nm, 1 < j < Nn, I < m,n < M + 1. 



(1) ^(2) (M+1) 



Theorem 2. Let T^v = 2iV^/3 ^^(^^^ ^ 2A^2/3 _ ^2/4 5 ^ R. For any M > 1, 

any sequence {A^m}m=i of positive integers, and any strictly increasing sequence {sm}m=i of 
nonpositive numbers with Sm+i = 0, 



, . . . , Sm+1,^m+1 
lim (Tat + Si, 6'aj^(si)^f • • • ; Tat, 6'aiv(sA/+i)^M+l' 



N 

Tdet Q -^''^'^^ 



where 9u{xi} = {xi + u}. 



Remark 2. This theorem may define an infinite particle system, in which any type of space- 
time correlation function is given by the above quaternion determinant. This quaternion 
determinantal system is the same as that derived in Forrester, Nagao and Honner 0, and it 
is inhomogeneous both in space and time. The spatial inhomogeneity is attributed to the fact 
that this system expresses the edge property of the infinite non-colliding Brownian particles. 
Thus, if we take the bulk limit, x\ — >• —00 with the position differences xf^^ — x^""* fixed, 
then the system should recover spatial homogeneity. It is confirmed by observing that the 
quaternion determinantal system given in Theorem 1 can be derived as the bulk limit of the 
system of Theorem 2, if we use the asymptotic expansion of the Airy function ()2.fij) 

~ ^1/2^^1/4 cos - as x ^ 00. 

On the other hand, keeping the spatial inhomogeneity, one can consider the limit Sm —00 
with the time difference Sn — Sm fixed, 1 < m,n < M. In this limit, V{sm,x] Sn,y) — > 0, 
T{sm,x; Sn,y) 0, and 

POO 

S{sm, X- y)^ dX e^'^-'-^'^/^Aiix + A)Ai(|/ + A). 
Jo 



Hence the off-diagonal elements vanish in the 2x2 matrix expressions of quaternion q™'"(x, y) 
and 

a'"'"(x,?/) 
a"'™(i/,a;) 



C{q'^'^{x,y)) 
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for 1 < m,n < M, where 

a"''"(a;,?/) = a(s^,x; s„,?/) 

dX e^'"-'"^^^/^Ai{x + X)Ai{y + A) if m > n 

d\ e^''""'''"^^/^Ai(a; + \)Ai{y + A) if m < n. 



Then the quaternion determinant TdetQ ^xj^^\ x^^, . . . , x^^^ j is reduced to an ordinary 

determinant det ^ ^xj^^, x^^, . . . , x^^^^ j with the elements a™'" ^a;-™''', x^"'' j . In this way, 

we will obtain the infinite particle system, which is temporally homogeneous but spatially 
inhomogeneous with the multitime correlation functions 

In particular, if we set iVi = A^2 = ■ ■ ■ = Nm = 1, then 

p'(si,{x«};...,SM,{xW}) = det (a'"'"(x('"),x("))). 

l<m,n<M 



This is the same as the system called the Airy process by Prahofer and Spohn in |2E]- (See 
also HI].) 

3 Quaternion determinantal expressions of the corre- 
lations 

In this section we give quaternion determinantal expressions for the correlation functions 
defined in fl2.5p along the procedure in [20] ■ From now on we consider the case N is even, for 
simplicity of notations. See 20^, for necessary modifications for odd case. For 1 < m,n < 
M + 1, define 

(3.1) = r r 7-'~(--;> ^-'"-(-•-) 

J~oo J^oo PT-t„[y,Z) Pt-uAV^W) 

where po{y,x)dy = 5x{dy). We introduce an antisymmetric inner products 

(/,^)™= I dx f dy F"^'"'{x,y)f{x)g{y), 
Jr. Jr 

and 

if, 9)= dx dy F^'\x,y)pt,{0,x)pt,{0,y)f{x)g{y). 
Jn Jn 

For /c = 0, 1, ... we consider the polynomials in x of degree k defined by 

(3.2) Rk{x) = z^" J2 ^k,H, zl 

j=i V^i/ 



where c, = = ^f^, 

(3.3) 



2 '^c^dkj, 



if A; is even, 

2"''cf {^j - 2(A; - l)4-2i}, if k is odd, 



and Hj{x) are the Hermite polynomials. They are monic and satisfy the skew orthogonal 
relations: 

{R2j,R2£+l) = —{R2t+li R2j) = ^j^jii 

(i?2j, R2i) = {R2j+i, R2e+i) —0, j, £ = 0, 1, 2, • • • , 



where 



r(j + |)r(j + i) ftlY'^'/' 



For m = 1, 2, . . . , M + 1, and A; = 0, 1, . . ., put 

(3.4) Rt^\x)= [ dy Rk{y)ptAO:y)Pt^-tAy:^)- 



Then we can prove the skew orthogonal relations 



2£+H -'T'2j /m — ^j^jt-, 

— /R^™) \ — n T — n 1 9 

\-n,2j , /m — \-'^2j+l' -'^2f+l/™ ~ J) ~ ^) 

for any m = 1, 2, . . . , M + 1. For m = 1, 2, . . . , M + 1, define 



(3.5) 



dy Rt\y)F"^^"^{y,x), fc = 0,l,2,... 



R 



Now we introduce the functions on R^, and S""*'", 1 < m,n < M + 1, given by 

(7V/2)-l ^ 



(3.6) 
(3.7) 
(3.8) 



I k 



k=0 

(iV/2)-l 



(2;,y) = 



fe=0 
{7V/2)-l 



fe=0 



Further we define 

(3.9) 
(3.10) 



5"'"(x,i/) = S^'-{x,y)-p,„.,Jx,y)l{m < n), 
i'^'''{x,y) = I"''^{x,y) + F"^'^{x,y). 
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Define the quaternions q"^'^{x,y), 1 < 171,11 < M -\- l,x,y E R so that these 2x2 matrix 
expressions C{q"^'^{x,y)) are given by 



C(g™'"(a;,y)) 



1 < ^ < AT^, 



Let M > 1 and {Nm}m=i be a sequence of positive integers less than or equal to A^. 
For xS7^ G R^, 1 < m < M + 1, we denote by Q (^xj^J, xjJJ, . . . , xj^^+^^^j the self-dual 

X]m=i X Sm=i quatemiou matrix whose elements are g™''" , x^"'' 
1 ^ J ^ 1 < wi, < M + 1. Then we show the following relation. 

Theorem 3. The multitime correlation function fl2.5p is written as 

(1) ,J2) (M+1) 



pI (ti, e^; . . . ; tM+i, ^M+t') = Tdetg (xSi;, X 



In order to prove the theorem, first we introduce the Pfaffian. For an integer and an 
antisjTiimetric 2N x 2N matrix A = (aij), the Pfaffian is defined as 



Fi{A) = Ffi<i<j<2N{aij) 
1 

~ m 



rry ^ sgn(cr)a^(i)^(2)a^(3)<j{4) ■ ■ ■ aa{2N-i)a{2N), 



where the summation is extended over all permutations a of (1,2,..., 2A^) with restriction 
a{2k — 1) < cr(2/c), A; = 1, 2, . . . , A^. If Q is an iV x iV self-dual quaternion matrix, then 

(3.11) TdetQ = Pf(^JC(Q)), 

where J is an 2A^ x 2A^ antisymmetric matrix with only non-zero elements 

J2k+l,2k+2 = —J2k+2,2k+l = 1, /c = 0, 1, 2, . . . , — 1. 

See, for instance, Mehta jTJj. 

For a function -i/;™-'" defined on R^ we denote the A^ x A^-matrices whose (i,j)-entry 

is ■?/;'^'" , x^"'' j by ip'^''^ ^xj^"*, xj^'' j , or simply by t/;"^'" for short. And we denote by 
^(m) j^^Mj the A^ X A^ matrix with (xJT^) = R^^\{xi), and by (^xj^^) that 
with i'^^Ni = ^'f^i{xi). Let L be the N x N diagonal matrix with Lj^j = ^r[(j_i) 
i = 1, 2, . . . , AT, and (^xj^^^) = L-ii?^'") (^xj^^^) . Then we have 



i-l)/2], 



(3.12) i?^'") (xJi^) Ji?(") (xSi^)^ = D"''^ (xS^U^") 
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As the first step of the proof of Theorem 3. We show that the multitime probabihty 
density defined in ()2.3|) is written as 



(3.13) pI (ti, ; . . . ; ^m+i, ^m+i) = TdetQ (x«, 



N 1 ■ ■ ■ 1 



(M+1) 



For simphcity of notation, here we give the proof of ()3.13|) for M = 2. It is straightforward 
to prove ()3.13|) for general M. Since 



sgn (jiN (^Sv^)) = Pfi<i<j<Af (sgn (^xf' - xf^ 



and sgn{y — x) = F'^' (x, y), we have 



(3.14) 



sgn ( /lAT f X 



(3) 



Pf [F^'^] . 



Noting that Rk{x) is the monic polynomial of degree k, we have 



hN (x(^)) = det 



X 



(1) 



i-1 



det ( Ri-i ( X 



„(i) 



and so 
(3.15) 



N 



l[pt, (0,xfM/ijv (x(i)) =det 



(1) 



R^'> X 



i=l 



Since det L = Uk=o~^ n = C{N,T,ti)-\ from and 



N 



(3.16) 



i=l 



C{N,T,h)l[pt, (0,a;fM (x«) = det (x^ 



I X^ , Ji.jy 



= Pf 
= Pf 

Then from ((231), (j3Ti|) and (j3l6|l we have 



Pf[D^.^]Pf[F3.3] n ^^et^ 



^« (x«) J^(^) (x« 



(-l)3^/2pf 



m=l 

Di'i o 



O -F3'3 

By basic properties of the Pfaffians, we have 



npf 

m=l 



o 



o 



Pf 



O 



npf 



-' m=l 



11 



Pf 



Pf 



Pf 
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~{Pts-t2, 
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Since G R< , ( x]y ) 7^ 0, and so det 
matrices 



.(1) 



7^ by ()3.15|) . Hence we can define 



.(1) 



-1 



which satisfies 



By repeating elementary operations, we see that the last Pfaffian equals to 



Pf 



r Di'i 




/)1,2 


(^2,1 


)T 


£)1,3 


(53,1)T 








_gl,2 


_/l,2 


_^1,3 


_jl,3 




(^1,2)T 


D2,2 






D2,3 


^3.2 






./2,1 


_g2,2 


_/2,2 


_^2,3 


_/2,3 






^3,2 




r 


^3,3 


(^3,3)T 


_53,1 




.J3,l 


_53,2 


_/3,2 




_/3,3 






- ^1,1 


^1,2 


^1,3 - 








_l)37V/2pf 


^2,1 


^2,2 


^2,3 












^3,1 


^3,2 


^3,3 
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where each A"^'"- = (A^'^) is a 2N x 2N matrix which consists of 2 x 2 blocks 



4Tn,n 



ij ji 



Qm,n j'tri. 
I m,Ti\ 



We can see that the above matrix A = {A^^' ) satisfies the relation A = JC{Q). Therefore, 
(jTT^ is derived from (jSZH). 

For square integrable functions and ip defined on R^, put (t)*ip{x, y) = z)tp{z, y)dz. 

Then we have 

^m.P ^ gP,m _ jm,p ^ jjp,n _ jjm,p ^ pp,n _ gm,p 

j-^m,p ^p,n j~^m,n (^m,p ^ jp,n ^m,p ^ pp,n jm,n 

S""'" * Pt„-t, = 5'"'", D^'P * pt^_t^ = D'"'", if p < n. 
Hence by simple calculation we see that 

q"''"'{z,z)dz = N, 
q^'^ix, z)q^'-{z, y)dz = g'"'"(a;, y) + ^'"'"(x, yHn, p) - k{p, m)g'"''^(x, y), 
where n{n,p) is a quaternion with 

1 - l(p < n) 

-l{n<p) 



R 



C{K{n,p)) 



Then by slight modification of Theorem 6 in j22] we have the following integral formula for 
any 1 < iV^ < A^, m = 1, 2, . . . , M + 1, 

1 uei^y 1 x^^ , . . . , x^^ , . . . , x^^^^^ I ax^^ 

R ^ ^ 

= (iv - iv„ + i)Tdetg (x«, . . . , x£Li, . . . , x^;^;;^)) , 

which is the generalization of the formula ()1.2p given in Introduction of the present paper. 
Successive application of the above relation yields Theorem 3. 



4 Expansion using Hermite polynomials 

In this section we show expansions of functions Ptn-tm^ -^i"*^ ^.nd by using Hermite 
polynomials Hk. Put 



7^ , 7n = ■ 



T \/ + 
V t-n 
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and r*^") = — logz„. By simple calculation we have 

g-(i™/2T)(x7c„)2g( 



A/27r(t„ - tm) 
( 



X exp 



(z//c„)-e-(^'"'-*™')(x/c„) 

]_ _ g-2{T(")-r{'"))^ 

\ 

for l<m<n<M + l. Using Mehler's formula |2] 

exp (-^^^^) = e-^V^(l-^^) f2^H,ix)H,iy) 



we will have the following expansions using the Hermite polynomials. For 1 < m < n < 
M+1, 



(4.1) 



2^e-|{l+7m)(x/c™)2g-|(l-7„)(j//c„)2 



VU^T - t 



><J2- — ^ — 

k=0 



hk 



X 



and for 1 < m < M + 1, 

Pti(0,a;)pt^_ti(x,?/) 



y/2^g-(x/ci)2g-i(l-7„)(i//c™)2 



iy/27rtitm(2T - ti) 

°° g-fe(T('")-T(l)) 

X XI ^ 

fc=0 



Then from ()3.2p . ()3.4p and the orthogonal relation of the Hermitian polynomials, we obtain 
(4.2) (x) = ^ Y.-^,,e-^ ( f 



j=0 



e 2 



|(l+7m){3;Vc„)2g-l(l+7,0{j//c„)2 



From the definition ()3.1|) and the expansion ()4.1|) we can obtain 
(4.3) F^'-{x,y) 

oo « 



V(2T-(,„)(2r-(, 



°0 °0 gfcT('")g£T(") 



i^^A. — i// - O/i, 



1/ 



where (■,■)* is the antisymmetric inner product defined by 



El 



oo •/ — oo 
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f{z)g{w) - f{w)g{z) 



Put Rl{x) = Yl^=o'^kjHj (^■^:^^^ . Then {Rl{x)} satisfy the following skew orthogonal 
relations 

(4.4) (i?2j, -RWl)* = ~{R*2e+lJ R*2j)* = ^*j^3h 

{r;^, r;,). = 0, r;,^,). = o, for j, £ = o, i, 2, . . . , 

where r| = 4h2iT{ci/2)^^'+\ We put 

{2^c{^5jki if ^ is even, 

2^(^)'{ci(^)'} \ if A;, J are odd and > J, 
0, otherwise, 



for nonnegative integers k and j. Then X]j=s /^^i'^i* ~ ^ks-, if < s < /c, and 

(4.5) =X:/?..^K^)- 

From the definition ()H.5|1 and the equations ()4.2|1 and ()4.Hp we have 

m j 



X 



i:i:v*(l:)-'(^'b)'^'vvT 

V27rT(2T - «„) £^ ft< VcJ 

Using the skew orthogonal relations ()4.4|1 . we show that for A; = 0, 1, 2, . . . 



. . p-^ 1+7™ Wcm V* ,n ^ p^^'"' / r \ 

(4.6) $r^(x) = , ^ V 

v/27rT(2T - 1„) ^J^^^ /^^ Vcn^y 

Using above expansions we show the following lemma. 
Lemma 4. For 1 < m, n < M + 1, 



(4.8) F"^''^(x,y) = ^ ^ [<^>^^(a:)<f£Vi(^) - <^2rii(^)<^'2:^(^/) 

fc=0 ^ 
00 ^ 

(4.9) r"'"(x,y) = E r Kr^(^)*2ti(^) - *2rli(^)*2? (1/) 



fc=Af/2 
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Proof. By (gTj) and the relation 

(4.10) rk 



1 / + \ 2fc+l/2 



27rT V2T-ti 

we have 



oo 



fc=0 ^ 

e-|{l+7m)(x/c,„)2g-i(l+7„){y/c„)2 oo 



v/(2T-tJ(2T-tO 



fc=0 



{oo ^vt™) / \ °o 

j=2k ^^'"Z ^=2fc4 



— ; Pi2k+lHi — 

2k+l ^ ^ " 



By fj4.4|) and ()4.5p the right hand side of the above equation equals to 

g-|(l+7m)(x/c™)2g-|(l+7n)(2^/c„)2 OO OO 



oo oo v(m) «^{n) / \ / „ \ 



where we have used ()4.H|1 . From the definitions ()H.7|1 and (jH.lOj) . ()4.9j] is derived from ()4.8j) . 

5 Proof of Theorems 

The following formulae are known for (jl.lj) [21 EH] • For m G R, 

(5.1) (^) 

(5.2) lu„(-l)¥/V.« (^) 

(5.3) hm 2-i<Av'< ^) = Ai(-«) 

Here we give the proof of Theorem 2 by using ()5.3|) . The proof of Theorem 1 will be easier 
and given by the similar argument using ()5.1|1 and ()5.2j) . 



u \ 1 

COSM, 
U \ 1 

sin u, 

'TX 
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Let = ^2T - exp |l/27„(x/c„)2 - iVr('")| and (""(x) be the quaternion with 

b"Hx) 



C(r(x)) 



1/6™ (x) 



For xj^"* G R^, l<m<M + l, we consider the transformation of the quaternions 
g,m,n x^") j g™'" ^X-'"^ xj""* j defined by 

We denote by Q (^xj^^\ xj^^, . . . , ^^NmS) ^'^^ self-dual N^. x Z]m=i quaternion ma- 

trix whose elements are g™'" ^Xj-™\ x^"'* j , 1 < < A^mj 1 < J < ^n, 1 < "m, n < M + 1. By 
the definition of quaternion determinants, the following invariance is established: 

Tdet O ("x^') x^'^ x^*^"^') ^ - Tdet O Tx^'^ x^™^ x^'^^'^ ^ 
1 uei i,^ ^x^^ , x^2 , . . . , y^j^^^^^^ J — 1 uei i,^ ^x^^ , x^^ , . . . , -X-atjv^+i J ■ 

Hence to prove Theorem 2 it is enough to show the following lemma. 

Lemma 5. Let T/v = 2N^^'^ and t„i = Tn + Sm, l<m, n<M+l. Then for any x, ?/ G R, 

(5.4) lim \ D"^'''{aN{sm) + x,aN{sn) +y) = 'D{s^,x;Sn,y), 

(5.5) lim 6'"(x)6"(?/)/"''"(a^(s^) + x, a7v(sn) + 2/) = i{sm, x; y), 

(5.6) lim \-^S"'''^{aN{Srn) +X,aN{Sn) +y) = S{Srn,X]Sn,y). 

N^oo [y) 

We start to prove this lemma by showing 

b'^(x) 

(5.7) lim ——pt^^t^{aN{Srn)+X,aN{Sn)+y) = V{Srn,X]Sn,y). 

By dHH) and the fact 

J-N 



for large A^, we have 



6™(x) 

-^;7^Pt„-t™(a7v(s^m) + x, aiv(s„) + y) 



•L- t e^'--'Ai (. + ^) Ai (, + ^) . 



p=— oo 
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where we have used fl5.3|) . Then we have ()5.7|) . 
From (jnSl), (Q, (ESI) and (glTj), we have 



with 



S-'-{x,y) = Sr{x,y) + ST''{x,y), 



cm,n/ \ ^"(y) (Ar_£)(T-(")-r('"') f \ f y \ 



oo 



^ BCiV/2-l-i) ,„_2j_,,,(„) / a: 

X 1. BiN/2-lf ^'-^ i^:; 



k=N/2 



where B(k) = ,^ ^' Since 

by the same argument to show ()5.7|) we have 

b'^(x) 

(5.8) hm ———S"''''{aN{srn)+x,aN{sn) + y) = S{s^,x-,Sn,y). 

()5.6|) is derived from ()5.7|1 and ()5.8|1 . 

From ()4.fij) and ()4.7|1 . by calculations with ()4.1()|1 . we have 

y»(x)r(»)4.K2,(i^)*5;l2p+i(») 



X 

A;=p 

From ()4.9|) we obtain ()5.5p by the same procedure as above. 
From ()4.2p . by calculations with ()3.3p and 



we have 



rkh'^{x)h'^{y) 



'1 ,(7V-2fc)r(™) \ (7V-2fc+l)r(") 
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y_ 



Using the fact that 

(jV-p)i/i2 f a„(s„) + y \ d 

for large A^, we obtain ()5.4p . This completes the proof of Lemma 5. 
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